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We numerically study quantum phase transitions and dynamical properties in the one-dimensional cluster model with 
several interactions by using the time-evolving block decimation method for infinite systems and the exact diagonaliza- 
tion. First, boundaries among several quantum phases of the model are determined from energy gap and each phase is 
characterized by order parameters and the entanglement spectrum (ES). We confirm that in the model with open bound¬ 
ary condition the degeneracy of the lowest levels in the ES corresponds to that of the ground states. Then, using the 
time-dependent Bogoliubov transformation with open boundary condition, we investigate dynamical properties during 
an interaction sweep through the critical point which separates two topological phases involving four-fold degeneracy in 
the ground state. After a slow sweep across the critical point, we observe spatially periodic structures in the string corre¬ 
lation functions and the entanglement entropy. It is shown that the periodicities stem from the Bogoliubov quasiparticles 
generated near the critical point. 


Introduction- The exploration of topological phases such 
as quantum Hall states, topological insulators, and spin liq¬ 
uid phases has been carried out intensively for three decades. 
Topological phases are characterized not by any local order 
parameters but by nonlocal order parameters or their emergent 
edge excitations. 1_8) Recently, it was found that entanglement 
can characterize topological phases. 9-12 ) Li and Haldane pro¬ 
posed the concept of entanglement spectrum (ES), which is 
obtained from the eigenvalues of the reduced density matrix 
of a subsystem. 13 ) They found that the ES of fractional quan¬ 
tum Hall states is similar to the energy spectrum of the low- 
lying excitations appearing at the edge of the system. Since 
then, the ES has been widely used as a tool to study topologi¬ 
cal phases. 14-23) 

One of the simplest models which show a variety of topo¬ 
logical phases would be the cluster model. 24-26 ) The ground 
state of the model is called the cluster state, 27 ) which is char¬ 
acterized by the non-local string order parameter. 26 ) In the 
case of open boundary condition, two zero-energy modes 
localized at each end of the system exist and four-fold- 
degenerate ground states result. The model is interesting also 
from the quantum-information perspective. In fact, one-way 
quantum computation and measurement-based quantum com¬ 
putation using the cluster state were proposed. 27-29 ) In addi¬ 
tion, the one-dimensional cluster model is expected to be re¬ 
alized in experiments of cold atoms on a zigzag ladder by 
introducing three-spin exchange interaction. 30 ) 

Dynamical properties of systems associated with phase 
transitions have been extensively studied for a long time. 31-33 ) 
Dynamics during a parameter sweep across the critical point 
especially has been investigated. 34-39 ) In phases characterized 
by conventional local order parameters, a universal relation 
called the Kibble-Zurek scaling 32,33) is known to hold for the 
dynamics of defect density. However, dynamics in topologi¬ 


cal systems depends on their topological properties and differs 
from the Kibble-Zurek physics, as pointed out in Refs. 4CM-2. 

In this letter, we focus on topological quantum phase tran¬ 
sitions and dynamical properties under an interaction sweep 
through a critical point associated with a topological phase 
transition. Specifically, we generalize the cluster model in one 
dimension by adding interaction terms with topological na¬ 
ture, which give rise to a quantum phase transition between 
two different cluster states. Our model would be interesting 
in the light both of quantum-information science and of com¬ 
peting topological phases. We map out the ground-state phase 
diagram of our model by calculating energy gap and char¬ 
acterize the phases by order parameters and the ES. In addi¬ 
tion, the dynamics during an interaction sweep across the crit¬ 
ical point by calculating energy gap is investigated. For slow 
sweep speeds, we observe, after passing the critical point, pe¬ 
riodic structures in the length dependence of string correlation 
functions and the entanglement entropy (EE). This breakdown 
of adiabaticity is due to the topological degeneracy of the ini¬ 
tial cluster state. We conclude that the periodicity is attributed 
to the Bogoliubov quasiparticles (bogolons). 

Model- The one-dimensional cluster model is defined by 

Hc = ~t jXZX ^ M <2’ (D 

i= 1 

where N is the system size and cr ? (a = x, y, z ) are the Pauli 
matrices at site Z. 24-26 ) With open boundary condition, we sup¬ 
pose <r^ +1 = cr^ +2 = 0 (a = x,y, z). The three-site interaction 
in Eq. (1) is called the cluster interaction or the cluster stabi¬ 
lizer in quantum-information science. In the ground state of 
the model dubbed the cluster state, the string order parameter 
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is called the string correlation function of length L. 1,3,26) A 
phase characterized by the non-vanishing string order param¬ 
eter is generally called the cluster (C) phase. Throughout this 
letter, we consider a generalization of the cluster model in one 
dimension to clarify the topological properties. The general¬ 
ized cluster model is defined by adding the Ising interaction 
<tJ<tJ +1 and another cluster interaction c r fc r f + 1 c r f +2 t0 the orig¬ 
inal cluster model: 

Hoc = Yj^- jXZX ^<i^ 2 +jYY ^M +jYZY ^<yi + 2y 

(3) 

Note that our model at J YZY = 0 is the same as the clus¬ 
ter Ising model studied in Ref. 26. When J YZY is dominant, 
there appears another topological phase, which we call the 
dual cluster (C*) phase. The phase is characterized by the dual 
string order parameter Oyzy = lim^^oo Oyzy(L ), where 

(L -2 \ 


Oyzy(L) 


°>2 


\i =3 


1 <-y L 
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is called the dual string correlation function of length L. 

The model (3) can be solved by using spinless fermion rep- 
resentation. 43) The original spin model (3) is transformed into 
a quadratic Hamiltonian: 


N r i 

h =Yj c l A ‘jCj + - {c]Bijc) + CiBjiCj) 


(5) 


hj= 

by the Jordan-Wigner transformation: 

i -1 i -1 

Ci= c i =n(-5K. ^ 

y=l 7=1 

where A is a real symmetric matrix and B is a real antisym¬ 
metric matrix with A^+i = J YY , A iti+ 2 = J xzx - J YZY , B i i+ \ - 
- J YY , B i i+ 2 = J xzx + J YZY . In general, quadratic Hamiltoni¬ 
ans (5) can be diagonalized by the Bogoliubov transformation 
as 


h= Yj e ^ ^- 5 )’ E *-°- (7) 

The ground state is given by the Bogoliubov vacuum |vac) 
satisfying r 7 ^ |vac) = 0 for all fi. 

The topological nature of the model (3) is clearly seen in 
the Majorana representation. 4 ^ The Majorana fermions {q} 
are defined by real and imaginary parts of the two fermion 
operators {cj and {cj}: 

c' 2,-1 = cj + Ci, C 2 i = i (■ Ci - cj), i= 1,2,..., N. ( 8 ) 
The standard anticommutation relations of {c,} and {cj } trans- 


Ci = cj, {cj, cj) = 28ij. (9) 

Using the Majorana fermions, the model (3) is rewritten as 

H GC = )fMc, (10) 

where c - (ci, c 2 ,..., c 2 aO T and M is a real antisymmetric 
matrix with M 2i -_ u? - +2 = - J YY , M 2/ _ u/+4 = - J YZY , M 2i , 2 i +3 = 
j xzx y/e can easily read off the number of ground states and 
that of zero modes from the above Hamiltonian. Suppose the 
model with open boundary condition. When either J xzx or 
J YZY is dominant, two zero-energy modes localized at each 
end of the system exist and result in the four-fold degeneracy 
in the ground state, as shown in Fig. 1 (a), (c). In contrast, 
when J YY is dominant, only one zero-energy mode exists and 
the ground state is two-fold degenerate, as shown in Fig. 1 (b). 
Thus, the number of zero-modes appearing at the edges of the 
system characterizes each phase. 46) 

Phase diagram- In Fig. 1 (d), we show the ground-state 
phase diagram of the model (3) and characterize each phase 
by order parameters and the ES by using the method given 
in Ref. 44. Here we impose the periodic boundary condition. 
By performing the Bogoliubov transformation, the model is 
expressed in the momentum space as 

H = Ak + kl'-tjl-k)’ A k = 2 ^4 + 6 l’ J ) 

0 <k<7T 

where > 0 is the excitation energy at the wave number k 
and 

€k = ( J xzx - J YZY ) cos 2k + J YY cos k , ( 12 ) 

6k = ( J xzx + J YZY ) sin 2k - J YY sin k. (13) 

In the following, we use J xzx as the energy unit and assume 
that J xzx is positive. The phase boundaries depicted by the 
thick solid curves in Fig. 1 (d) are determined by the condition 
A^ = 0 at a certain k. 

In order to characterize each phase separated by the curves, 
we calculated the order parameters with the time-evolving 
block decimation method for infinite systems (iTEBD). 47,48) 
This method is advantageous over the exact diagonalization 
for the computation of the order parameters with no bound¬ 
ary effect. Figure 1 (e) displays the string order parameter (2) 
and ferromagnetic order parameters O aa = lim^co ^ cr®cr 
(<a = x,y) along the thick dotted line in Fig. 1 (d). Note that 
we practically take L = 200, at which the order parameters 
safely converge to their thermodynamic values. We first con¬ 
sider J YZY = 0. When \J YY / J xzx \ becomes large, conventional 
ferromagnetic (F^) and antiferromagnetic (AF (},) ) orders in 
the y-direction appear, which are respectively characterized 
by the ferromagnetic and antiferromagnetic order parameters 
defined by O yy and lim L ^oo(-l) L (cr^(T^). When \J YY /J xzx \ 
becomes relatively small, on the other hand, the C phase char¬ 
acterized by the string order parameter (2) appears. Next, var¬ 
ious phases appears due to non-zero J YZY are found. For pos- 
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Fig. 1. (Color Online) (a)-(c) Schematic representation of the interactions 
in Eq. (3) by the Majorana language, (a), (b), and (c) depict the first, second, 
and third terms of the Hamiltonian (3). Non-interacting Majorana fermions 
enclosed in dotted line appear at the ends of the system, (d) Phase diagram of 
the generalized cluster model (3) for J xzx > 0. On the thick solid curves, the 
excitation gap vanishes at a certain k. C, C*, F, and AF represent cluster, 
dual cluster, ferromagnetic, and antiferromagnetic phases, respectively. The P 
phase cannot be characterized by string and (anti) ferromagnetic order param¬ 
eters. The superscript represents the direction of the order. Each phase is de¬ 
termined by order parameters calculated with the iTEBD. Along the thin dot¬ 
ted line and the dashed line, we calculated the ES in Fig. 2 (b), (c). (e) String 
order parameter Oxzx and ferromagnetic order parameters Oxx,Oyy in the 
x, y-directions along the thick dotted line in Fig. 1 (d) (J YZY /J xzx = 0.5) by 
using the iTEBD with bond dimension x - 60. 


itive J YZY , the AF and F phases ordered in the v-direction (la¬ 
beled respectively by AF (x) and F w ) appears. For negative 
J YZY , on the other hand, we find the paramagnetic (P) phase 
which cannot be characterized by the string and the (anti) fer¬ 
romagnetic order parameters. Similar behavior was reported 
in the cluster-XT model 49) and in the transverse Ising model 
in frustrated lattices. 50) 

Let us interpret the phase diagram obtained above in terms 
of entanglement. Here we suppose open boundary condition 
and use the exact diagonalization. We divide the entire system 
into two subsystems A and B symmetrically around the cen¬ 
ter with the length of A being L su b (Fig. 2 (a)). We calculate 
the eigenvalues {A v } of the reduced density matrix p A of A 
defined by tracing out the subsystem B: p A = Tr B p, where p 
is the density matrix of the ground state of the entire system. 
The ES is defined as = - In d y . 13) First, we focus on the thin 
dotted line depicted in Fig. 1 (d). As shown in Fig. 2 (b), the 
degree of degeneracy in the lowest value of the ES stays con¬ 
stant in each phase. We see that its number is four, one, and 
two for C, P, and AF phases, respectively. Similarly, focusing 
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Fig. 2. (Color Online) (a) Schematic of subsystems A and B. (b) ES with 
N = 503 and L sub = 249 for J YZY /J xzx = -0.5 indicated by the thin dotted 
line in Fig. 1 (d). The degeneracy in the lowest levels is four, one (no degen¬ 
eracy) and two in C, N, and AF phases, respectively, (c) ES with N = 503 
and L sub = 249 for J YY /J xzx = 0 indicated by the dashed line in Fig. 1 (d). 
The degeneracy in the lowest levels is four in both C and C* phases. 


on the dashed line depicted in Fig. 1 (d) the degree of degen¬ 
eracy is four in both C and dual C* phases, as shown in Fig. 2 
(c). At the critical points, the degeneracy structure disappears. 
We also calculate the topological invariant, called the winding 
number, 5 ^ by using the method explained in Ref. 52. We con¬ 
firm that the winding numbers of C, C*, AF, F, and P phases 
are two, two, one, one, and zero, respectively, which corre¬ 
spond to the number of Majorana fermions appearing at the 
edge of the system (not shown). From all of these results, the 
degree of degeneracy is equal to the number of the Majorana 
zero-modes existing at the ends of the system. We thus corrob¬ 
orate that for the generalized cluster model (3) the ES reflects 
the fictitious edge modes appearing at the position where we 
cut the system as in the case of other topological phases stud¬ 
ied in previous works. 13-23) 

Critical sweep- Now we turn to the dynamics during an 
interaction sweep across the critical point between C and C* 
phases with open boundary condition. Let us consider the fol¬ 
lowing time-dependent Hamiltonian: 

H(t) = -J xzx 2 cr^ +1 of +2 + J(t) LX°1 + A+2’ (!4) 

i= 1 1=1 

where the interaction parameter changes linearly during the 
sweep time r as 

J(t)/J xzx = 2th, o < t < t. (15) 

Here, the eigenenergies in Eq. (7) are labelled in ascend¬ 
ing order; E\ < E 2 < • • • E N . Since E\ = £2 = 0 ex¬ 
cept J xzx = J YZY , tj[ |vac), rj 2 |vac), and rj^rjl |vac> are the 
ground states as well as the vacuum state |vac). The initial 
state is prepared in the Bogoliubov vacuum corresponding the 
ground state of the Hamiltonian with J(t) = 0. We calculated 


3 




































































the length dependence of the dual string correlation function 
OyzyW) and the EE S (€) by using the time-dependent Bogoli- 
ubov transformation, 53,54) where € is defined as Fig. 3 (a). We 
define the EE as the von Neumann entropy of p{l) which is 
the reduced density matrix of the subsystem with the length 
that is, S (£) = - Tr p(€) lnp(T). Figure 3 (b) shows the length 
dependence of the dual string correlation function and the EE 
in the final state (t = r) with r = 25, 50, 100, and 200. For 
large r, a quadruple-periodic structure in the length depen¬ 
dence is observed in both quantities. We calculated the ex¬ 
pectation values of the number of bogolons in the final state. 
Figure 3 (c) indicates that only third bogolon is dominant for 
larger r. At t = r/2, the instantaneous Hamiltonian represents 
the critical system. Since the energy levels of low-lying ex¬ 
cited states approach the ground states near the critical point 
as in Fig. 3 (d), transitions from the ground state to the excited 
states would occur for finite r. From this fact, one may specu¬ 
late that the similar periodic structures in the dual string order 
parameter and the EE come from the bogolons generated near 
the critical point. 

To substantiate this, we calculated the dual string correla¬ 
tion function for excited states. A similar periodic structure is 
actually observed in the length dependence of the string cor¬ 
relation function in the excited state with a single bogolon, 
as shown in Fig. 3 (e). This period corresponds to the wave 
length of the third bogolon. As mentioned before, the bogolon 
with zero energy corresponds to the ground-state degeneracy. 
Actually, as shown in Fig. 3 (e), the dual string correlation 
function in the bogolons with zero energy is the same as that 
in the Bogoliubov vacuum. 

Conclusion- We have studied the one-dimensional cluster 
model with several interactions to understand the dynamics in 
the topological phase. First, we have determined the ground- 
state topological phase diagram (Fig. 1) and characterized 
each phase by the order parameters. In particular, we have 
found a direct quantum phase transition between two differ¬ 
ent cluster phases. We have confirmed that the degeneracy in 
the lowest levels of the ES in each phase corresponds to the 
number of the Majorana fermions existing at the edges of the 
system (Fig. 2). Second, we have investigated the dynamics 
during the interaction sweep with finite speeds across the crit¬ 
ical point separating the two cluster states. The periodicity in 
the length dependence of the dual string correlation function 
and the EE has been observed (Fig. 3). We have reproduced 
similar periodic structure by using the excited states and veri¬ 
fied that the periodicity stems from the bogolons excited when 
the system is close to the critical point. It would be interesting 
to see the results from the viewpoint of topological block¬ 
ing. 41 ’ 42 ) 
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Fig. 3. (Color Online) (a) We calculate the dual string correlation function 
of length € and take € adjacent sites as the subsystem to calculate the EE. (b) 
The length dependence of the dual string correlation function and the EE of 
final state (t = r) with N = 101 and r = 25, 50, 100, and 200. As r increases, 
a quadruple-periodic structure is clearly observed, (c) The expectation value 
of the number of bogolons in the final state. This indicates that only third 
bogolon is dominant for larger r. (d) Energy spectrum of low-lying excited 
states. There is a level crossing at the critical point at J YZY /J xzx = 1. (e) The 
dual string correlation function in the states with the bogolon at J YZY /J xzx = 
2. We can see a quadruple-periodic structure in the length dependence for 
third bogolon excited state. The dual string correlation function in the state 
with the bogolon having zero energy is the same as that in the Bogoliubov 
vacuum. 
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